Chaotic Hamiltonian systems revisited: Survival probability 
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We consider the dynamical system described by the area-preserving standard mapping. It is 
known for this system that Pit), the normalized number of recurrences staying in some given domain 
of the phase space at time t (so-clled "survival probability") has the power-law asymptotics, P(t) ~ 
t~ v . We present new semi-phenomenological arguments which enable us to map the dynamical 
system near the chaos border onto the effective "ultrametric diffusion" on the boundary of a tree- 
like space with hierarchically organized transition rates. In the frameworks of our approach we have 



estimated the exponent v as v 
rotation number. 



I. INTRODUCTION 



ln2/ln(l + r g ) 



1.44, where r„ 



(v5 — l)/2 is the critical 



In this paper we propose a new estimate for the 
" Poincare recurrences" (or "survival probability"), P(i), 
for the standard mapping in the vicinity of the chaos bor- 
der [l|, Q • Our consideration explicitly exploits the fact 
that the phase space of the standard mapping near the 
border separating chaotic and integrablc behaviors has 
self-similar scale-invariant structure consisting of hier- 
archical set of metastable islands of integrability ("can- 
tori") embedded in the "chaotic sea" Q. To be specific, 
we consider the area-preserving standard mapping 



Vt+i =Vt- K/ (27r) sin27ra; t 
x t +i = x t + yt+i mod 1 



(1) 



For K £ [0, K g [ the phase space of the system has 
disjoined islands of integrability which is destroyed as 



K — > K„ from below, where K„ 



0.97163540631 



Above the critical value K g the behavior of the system 
is less universal: many invariant Kolmogorov-Arnold- 
Moser (KAM) tori disappear, but still some islands 
of metastability survive around the biggest resonances. 
However, we should emphasize that the reorganization of 
the phase space above the value K g has no consequence 
for our consideration since we do not touch the region 
K > Kg and are interested in the survival probability 
only when K approaches K g from below. 

Our consideration of the survival probability is semi- 
phenomenological, that is we relay only on the measur- 
able "macroscopic" characteristics acquired in course of 
the iteration of the map ([T]). To be precise, we rely on 
the following well-established and confirmed facts: i) The 
number of principal resonances follows the Fibonacci se- 
quence when K /'Kg [3, [f|; ii) The generic behavior 
of phase trajectories is as follows: the phase trajectory 
stays in the vicinity of some resonance (low-flux can- 
tori) and then rapidly crosses the chaotic sea until an- 
other metastable low-flux cantori is reached [1, @; iii) 
The phase space of the standard mapping is self-similar 



being usually represented by a binary (i.e. 3-branching) 
Cayley tree [1, 0, S|; lv ) The survival probability has 
power-law asymptotic behavior,_P(i) ~ t~ v (for the first 
time this has been shown in 0, Hj] ) . 

Remind that survival probability is the normalized 
number of recurrences ((T|) which stay in some given do- 
main of the phase space at time t. Different research 
groups present different arguments for estimates of v, 
typically 1 < v < 3. The most intriguing contradiction 
concerns the discrepancy in reported values of v. The 
numerical simulations @, Q demonstrate v w 1.4 -r- 1.5, 
while almost all known analytic constructions give essen- 
tially larger exponents: v « 1.96 in Q and v w 3.05 
in [§[. The scaling analysis Q valid just near the chaos 
boundary gives v = 3. Th e sp ecial attention should be 
paid to the recent works |l0l-fl2l]. In [l(| the authors 
demonstrate that by an appropriate randomization of the 
" Markov tree model" proposed in B, B one can arrive at 
the value v w 1.57. The works [by, [l2j claim v = 3 for 
sticking of trajectories near K g in agreement with [f| and 
v = 3/2 for trapping of chaotic trajectories in the vicinity 
of cantori for K pa 2£/7r (where £ is a nonzero integer). 
The similar exponent, v = 3/2, was also obtained for a 
standard map in the work [1 31 ] . There is a point of view 
that the value v rj 1.4-5-1.5 corresponds to an intermedi- 
ate behavior of the system which is not yet reached the 
statio nary regime - see the corresponding discussion in 

The main aim of the present letter is to present some 
simple arguments in favor of the statement that the value 
v rj 1.4 -T- 1.5 could be the actual decay exponent of the 
survival probability P(t) for t — > oo without additional 
randomization of local transitional probabilities. Let us 
emphasize once more that in our approach we get rid 
of microscopic consideration of the detailed structure of 
quasiperiodic orbits and corresponding fluxes, but char- 
acterize the phase space of our system for a given cou- 
pling constant, K in the vicinity of the critical value K g , 
just by the hierarchy of resonances. 

According to 0, [13, [3 , the structure of the critical 
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KAM curve is determined by arithmetic properties of the 
rotation number, r, in the continued fraction representa- 
tion: 

r = —j = [miTO2TO3 ...] (2) 

mi H j 

m 2 H ■ 

m 3 + ... 

The ns best convergent to the rotation number r is r n = 
\m\rti2 ■■■ra n \ = p„/q n with mi = m2 = ... = m„ = 
1, where q n is the Fibonacci number. Recall that the 
Fibonacci numbers satisfy the recursion relation q n +\ = 

q n + q n -i; qo = o,qi = l. 

For K /• Kg the periodic trajectories with rotation 
numbers r n determine the structure of the phase space, 
converging to the critical boundary curve with r g (sec 
(3|). In the limit n — > oo (i.e. for q n — > oo) the phase 
space becomes self-similar with the scaling factor s n = 
q n /q n -i -> 1 + r g w 1.618, where r 9 = [111...] = (\/5 - 
l)/2 is the "golden mean". The convergents r n (n is 
fixed) characterize the positions of unstable fixed points 
of resonances for a given value of a coupling constant, K. 

According to [f|, the average local exit time, r n , from 
a given scale n can be estimated as t„ ~ \r g — r' n \ 2 /D n , 
where r' n is the ns convergent of the critical rotation num- 
ber, r g , and D n is the local diffusion coefficient. Since 
\r g — r' n \ ~ q~ 2 and D n ~ g ? 7 5 , one gets the estimate for 
the local exit time: 

T n ~ Sn (3) 

When n — > oo (i.e. when if — > 2£" ff ) one finds for the 
exit time, r, from the scale n the following asymptotic 
behavior: 

r n ~(l + r s )" = (1.618)" (4) 

To summarize, the following two facts Q constitute the 
basis of our semi-phenomenological consideration: 

• When approaching the chaos boundary from below, 
the average local exit time, r n , from the metastablc 
island of the hierarchy scale n is proportional to the 
number of periodic orbits, q n (see ([3])). This fact 
reflects the "local" structure of our phase space; 

• When n — > oo, the local exit time, r„, grows ex- 
ponentially with the hierarchy scale n (see (@J). 
This fact reflects the "global" structure of our phase 
space and allows to construct the estimate for the 
survival probability. 



II. DYNAMICS IN THE HIERARCHICAL 
TREE-LIKE PHASE SPACE: DIFFUSION ON 
THE BOUNDARY VS DIFFUSION IN THE BULK 

The hierarchical construction of self-similar sets im- 
plies that new smaller domains being properly magnified 



(rescaled) with the limiting scaling factor s g , coincide (in 
the statistical sense) with the former ("parent") domains. 
This construction refers implicitly to the tree-like geom- 
etry. The tree-like hierarchical geometry of the phase 
space of the standard mapping has been proposed in the 
seminal analytic works 0, 13 • The authors have supposed 
that the states of the system are the regions bounded by 
the low-flux cantori. Each state consists of an infinite 
hierarchy of low-flux cantori of smaller scale. Between 
any two adjacent cantori there are possible many other 
sub-hierarchies (or " island chains" in the terminology of 
Q). However only one such "island chain" in each hi- 
erarchical level is considered. Thus, the topology of the 
full phase space can be represented as a following dia- 
grammatic hierarchy, where the hierarchy depth (level) 
is labelled by the same index n, appeared already in J2]) 
as the cutoff in the continuous fraction expansion: 




Knowing the local transition rates between neighbor- 
ing states from microscopic computations, the authors 
in @, [|| have considered the random walk on the 3- 
branching tree and have derived the corresponding infi- 
nite hierarchy of recursive relations for the survival prob- 
ability. The plausible conjectures about the closure of 
this hierarchy and subsequent analysis of the dominant 
contribution to P(t) allowed to extract the decay expo- 
nent v « 1.96 in [7|. However this value is still rather 
far from numerically obtained exponent v ~ 1.4 -j- 1.5. 
Below we formulate an alternative point of view on dy- 
namics on hierarchical landscapes enabling to find the 
value of v much closer to known numerical value of the 
exponent v. 

The hierarchical description sketched above is fully 
consistent with the consideration of a "diffusion in a 
mountain landscape" [l9| appeared in a generic descrip- 
tion of "diffusion in hierarchies" regardless the subject 
of the chaotic Hamiltonian dynamics. For heuristic con- 
sideration, let two sequences of real numbers, A n and 
V n , be correspondingly the sizes of landscape valleys 
(basins) and the heights of passages (energy barriers) be- 
tween these valleys with respect to some reference (en- 
ergy) level. The basin sizes and the barrier heights can 
be introduced iterativcly. Suppose that some dynami- 
cal system is located in a basin Aq at the initial time 
moment. During the time t\ the system overpasses the 
lowest available height Vi and the probability to find the 
system in some state distributed initially on Aq, relaxes 
to the larger basin A\ D Aq. Inductively, if the prob- 
ability to find the system is located in the basin A n at 
time t n , then at time t n+ \ it relaxes into a larger domain 
A n+ i D A n by surmounting the lowest available height 
V n . Since by construction the n's basin hierarchically in- 
cludes all basins A„_i D A„_2 D ... D Aq, there is no 
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difference between the waiting time in the n's basin and 
the total time from the beginning of the dynamical pro- 
cess. Thus, the survival probability under diffusion on 
hierarchies depends essentially only on scaling of basin 
sizes, A n , and barrier heights, V n , and exactly this fact 
allows one to use the tree-like spaces with alternative 
(Archimedean or non- Archimedean) metrics (see, for ex- 
ample, 0, mi). 

Now we are in position to describe the main idea of the 
present work. In our description the state of the system 
for some value of K near the chaos border is uniquely 
characterized by the number of quasiperiodic orbits, q„ at 
the hierarchy level n. Thus, the states are parameterized 
by the Fibonacci numbers, q n . 

We consider the dynamics of the system as the transi- 
tions between different states for given value of n. This 
is the key difference with the former description of 0, H| 
schematically outlined above. Namely, in the former ap- 
proach the authors have considered the local random 
walk in the bulk of the 3-branching Cayley tree with the 
transitions between neighboring states belonging to two 
neighboring levels of hierarchy. Such a construction sug- 
gests an Archimedean metric on the tree-like space. To 
the contrary, we propose to consider an effective Markov 
dynamics on the boundary of the 3-branching Cayley tree 
truncated at the hierarchy level n (n^> 1). We allow for 
long distance jumps which appear with the probability 
prescribed by the limiting scaling factor s g only between 
the states of the same hierarchy level n. This construc- 
tion implicitly suggests non-Archimedean (ultrametric) 
space of states. 

To be precise, our Markov process takes place in an 
effective "energy landscape" constructed in the follow- 
ing way. We consider the phase space up to the scale n 
meaning that we regard the low-flux cantori (metastable 
islands) of scale n as local (possibly degenerated) min- 
ima of an energy landscape. Introduce now the selfsimi- 
lar scale-invariant structure of the basins of local minima 
hierarchically embedded into each others. Namely, each 
larger basin of minima consists of smaller basins, each 
of these consists of even smaller ones etc. Since the hi- 
erarchy level n is chosen arbitrary, each local minimum 
(i.e. low-flux cantorus of scale n) could be (and should 
be) understood as a basin as well, containing again the 
scale-invariant hierarchy (" subtree" ) of basins of smaller 
scale. Such a hierarchy does not conserve the structure 
of phase space as the islands of stability in the chaotic 
sea, but preserves the tree-like factorization of islands of 
smaller scale out of the islands of larger scale as K /*• K g . 
Note that in the common description of tree-like factor- 
ization of low-flux cantori, the value of n is counted from 
the root of the tree, meaning that larger and larger val- 
ues of n correspond to metastable islands of smaller and 
smaller scales. 

To specify our long-distance-jump Markov process in 
terms of the transition rates between the states (local 



minima), we construct the hierarchy of barriers between 
the basins of minima. Namely, larger basins are sepa- 
rated by higher barriers, while embedded smaller basins 
are separated by lower barriers. The Cayley tree can be 
regarded as a hierarchical " skeleton" of energy landscape 
(but not as a space of states as in 0,3): the bulk ver- 
tices of the tree parameterize the hierarchy of basins of 
local minima hierarchically embedded into each others; 
the same vertices parameterize the barriers separating 
the basins. 

The crucial requirement to our Markov process is as 
follows: the transition probability (per time unit) be- 
tween any two local minima is determined by the maxi- 
mal barrier separating these minima. This requirement is 
equivalent to the strong triangle inequality. Thus, basin 
sizes and barrier heights can be expressed in terms of 
ultrametric distances between local minima. 

In conventional description of ultrametric spaces ac- 
cepted in applications of p-adic mathematical analysis 
[22j to dynamics on energy landscapess, the ultramet- 
ric distances between local minima arc labeled in such 
a way that smaller transition rates (i.e. higher barriers) 
correspond to larger ultrametric distances. Therefore, 
while specifying ultrametric distances similar to the p- 
adic norm, p 7 , 7 = 1, 2, n (p is the prime number), we 
define the index 7 which numerates the hierarchy levels 
of Cayley tree in the direction opposite to the one of in- 
dex n. In this construction 7 = 1,2, ...n is counted from 
the boundary of the tree up to the tree origin, namely, 
there is a hierarchy of ultrametric distances between the 
states of scale n possessing the values p 1 ,p 2 , ...,p n . 

The hierarchical landscape with growing barriers corre- 
sponds to the fact that the average local exit time t„ from 
the states of scale n grows with n: larger the tree (i.e. 
deeper the hierarchy), longer the exit time averaged over 
all states on the tree boundary. This condition does not 
yet defines completely our energy landscape. We should 
satisfy another important requirement, viz, the average 
local exit time r„ from the states of scale n grows as 

r n ~ (q n /q n -lV * s n g (5) 

for n — > 00, where q n is the n's Fibonacci number and 
s g w 1.618. This means that the energy landscape, being 
expressed in terms of ultrametric distances between the 
states of scale n, meets a specific behavior of the average 
local exit time, r„, with n, and the survival probability 
as well, through the relation to the Fibonacci numbers. 
This is another key ingredient of our approach. 



III. THE TREE LIKE GEOMETRY AND THE 
FIBONACCI NUMBERS 

The procedure of construction of such a hierarchical 
landscape is described below. Anyway, to be specific, 
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in what follows we shall always keep in mind the 3- 
branching Cayley tree, C, as an example of an ultrametric 
space. 

The Fibonacci numbers have natural relation to the 
ultrametric geometry since they are connected to some 
discrete symmetries of the hyperbolic space. In order to 
substantiate our construction, it seems to be instructive 
to demonstrate briefly how the Fibonacci numbers ap- 
pear in ultrametric geometry and how they are connected 
to a 3-branching Cayley tree. Take the upper complex 
half-plane z = x + iy > and consider the zero-angled 
curvilinear triangle ABC bounded by two vertical lines 
AC, BC and a semi-circle AB leaned on the real axis 
as shown in the figure FigJTJt. Tessellate now the upper 
half-plane strip Imz >0;0<Rcz<lby the images 
of the triangle ABC . Two subsequent steps are shown 
in Fig[Tb,c. These images are obtained by sequential in- 
versions (fractional-linear transformations) of the initial 
triangle ABC. The images of all vertices of the triangle 
ABC lie on the real axis Im z = and the coordinates of 
the vertices of neighboring triangles satisfy the following 
composition rule 



X n +2 



Pn+2 




p n+ l dof Pn + Pn+1 



9n+l 



(6) 




A B'C'A' 
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1 3 2 3 1 

(c) 
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P„ + Pn» 
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In In+ln* 

(d) 



Pn 



Figure 1: Inversions of zero-angled triangles: a), b), c) - 
three subsequent stages of tessellation of the half-strip Im z > 
0; < Re z < 1 by the images of the triangle ABC; d) - 
composition rule for the coordinates of the vertices. 

Define now the 3-branching Cayley tree isometrically 
embedded in the upper half-plane strip Imz > 0; < 
Re z < 1 by connecting the centers of neighboring images 
of zero-angled triangles by arcs being parts of semicircles 
leaned against the real axis - see the FigJ5] Recall that 
the embedding of a Cayley tree C into the metric space 
is called " isometric" if C covers that space, preserving all 
angles and distances. For example, the rectangular lat- 
tice isometrically covers the Euclidean plane E{x, y} with 
the flat metric ds 2 E = dx 2 +dy 2 . In the same way the Cay- 
ley tree C isometrically covers the surface of the constant 
negative curvature, H. One of possible representations of 
%, known as a Poincare model, is the upper half-plane 
Im z > of the complex plane z = x + iy endowed with 
the metric ds% = (dx 2 + dy 2 )/y 2 of constant negative 
curvature. The composition rule ([5]) defines the coordi- 
nates of corresponding triangle vertices, {x n }. Hence, 



the rule (3) defines a set of rational numbers parameter- 
izing all bulk vertices of the Cayley tree of n hierarchical 
levels. The Fibonacci sequences appear as the subsets of 
{x n } corresponding to alternating left-right (or symmet- 
ric right-left) sequences of reflections of triangles (these 
sequences are marked in FigJ2]in boldface). For example, 
two "zigzags" starting from the point 0\ correspond to 
two "principal" Fibonacci sequences, f + {0\) = [1111...] 
and f~(0\) = 1 — [1111...]. The "secondary" zigzags, 
/(O2), /(O3) and /(O4) starting from the points O2, O3 
and 04 correspondingly, have the continued fraction ex- 
pansions: f(0 2 ) = 1 - [121111...], /(0 3 ) = [112111...] 
and /(O3) = [131111...]. Generically, all "zigzags" have 
the following continued fraction representation: 

{[1 7B2 777,3 •■• rrik 1111...] forl/2<a;<l 
1 - [lm 2 7713 ...m fe 1111...] for0<a;<l/2 

(7) 

with k first arbitrary numbers 7712, 7713, ...,m& (mi = 1) 
are followed by the "Fibonacci tail" of "1". 




1 1_ 2 1 3 2 3 
5 4 7 3 8 5 7 



1 lit 2 5 3 4 

2 7583 745 



Figure 2: Few subsequent inversions of the zero-angled trian- 
gle in the strip Imz > 0; < Re 2 < 1 is shown). Inversions, 
corresponding to Fibonacci sequences are shown in boldface. 

Washing out the metric structure of the set of 
fractional-linear transformations depicted in Fig(2] and 
leaving topological structure of the corresponding 
" zigzags" on the Cayley tree we arrive at the relation be- 
tween the Fibonacci numbers and an ultrametric space 
shown in Fig EH 

It should be emphasized that the ultrametric graph 
in FigJ3] designates schematically the hierarchically or- 
ganized transition rates (energy barriers) between the 
states on the tree boundary, while the bottom gray points 
labeled by rational numbers, p n /q n , parameterize the 
bulk vertices of the ultrametric tree, and hence, they cor- 
respond to the transition rates over the barriers. Again, 
the dynamics occurs only at this boundary, but not in 
the bulk of the tree. 



IV. SURVIVAL PROBABILITY 



Let us come back to the computation of the survival 
probability, P(t). Our main conjecture is as follows. Ex- 
tend the set of states of the dynamic system and suppose 
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Figure 3: Topological structure of the graph obtained by the 
successive applications of fractional-linear transformations in 
FigfJ] "Principal" and "secondary" Fibonacci sequences are 
marked in boldface (see the explanations in the text). 

that all fractions p n /q n parameterize some quasiperiodic 
orbits. The transition rates between any two different 
states M and N are defined by the relative "distance" 
between these states measured in number of successive 
reflections (as in Figj2|) necessary to superpose the points 
M and N. One might thought about our system in the 
following terms. Take a 3-branching Cayley tree up to 
some hierarchical level (generation) n shown in FigJ3] 
Parameterize all the bulk vertices of the tree by the sets 
of rational numbers p n /qn following the reflection de- 
scribed above, and consider the "ultramctric diffusion" 
on the boundary - a Markov process with the transition 
rates encoded in block-hierarchical kinetic matrix, known 
as "Parisi transition matrix" - see [2l|, (23|. Thus, our 
jump-like Markov process is defined by the Parisi kinetic 
matrix whose matrix elements are encoded by the set of 
rational numbers p n /q n - 

The survival probability, P(t), is the probability to 
find the system in the initial state after t jumps on the 
boundary of the Cayley tree (with hierarchically orga- 
nized transiti on p robabilities) . According to the works 
[2~i| (see also |25| for the transparent geometrical inter- 
pretation), the function P(t) consists of additive contri- 
butions from all possible directed paths on the p-adic 
tree (recall that in our case p = 2): 

P(t,T) = (p-1) P" 7 e A ^'+p- r (8) 

{7J} 

The indices 7 and j label correspondingly the hierarchical 
level of the tree (1 < 7 < T) and the specific point in the 
hierarchical level 7 (1 < j < p r ~ 7 ). Pay attention that 
now 7 is counted from the boundary of the tree towards 
the root point (see the discussion above) . The eigenvalues 
A 7J , are defined via the following construction 

r 

7' =7+1 



where the sum denoted by E runs along the tree from 
some vertex point labelled by the pair of indices (7, j) to 

(7) 

the root point 0±, and q-y is the transition probability 
corresponding to the state labelled by (7, j). 

The whole variety of directed sequences running 
from the hierarchical level 7 = 1 to the root point 
0\ is bounded by two "limiting" trajectories "loga- 
rithmic" and "linear". The "logarithmic" is o\ Q% = 
{1/5,1/4,1/3,1/2}, and the "linear", being the "prin- 
cipal" Fibonacci one, is cr Un = {3/8,2/5,1/3,1/2}. The 
denominators in the " logarithmic" sequence grow linearly 
with 7, g 7 = 7, while the denominators of the "linear" 
sequence grow exponentially, q-y ~ (1.618) 7 . The no- 
tations "logarithmic" and "linear" come from the fact 
that V(j,j) ~ — lng 7 can be considered as the effec- 
tive dimensionless local height of the potential barrier in 
the point (7, j). The "logarithmic" landscape is associ- 
ated with the "logarithmic" sequence, for which one has 
V(7) ~ In 7, while the "linear" landscape is associated 
with the "linear" sequence, for which one has V(j) ~ 7. 

Taking into account that: i) the eigenvalues © of the 
transition matrix arc given by the weighted sums along 
different directed paths on the Cayley tree, and ii) the 
survival probability does not depend on the multiplicity 
of paths on the tree with the same sequence of transi- 
tion probabilities, we can directly use the results of [24| 
for survival probabilities in logarithmic and linear land- 
scapes. 

The survival probability for logarithmic and linear 
landscapes reads (see (24| ) 

j e ~ t/ln2 for q-y = 2^ ^ (logarithmic) 
P ^ [C r 1/a for q y = 2~ 7 2~ Q7 (linear) 

(10) 

where c = T (£ + 1) (-r p=2 (-aV) 1/a and r p=2 (...) is 
the p = 2-adic T-function (see [24| for details). 

The appearance of the factor 2~ 7 in the transition 
probabilities g 7 in (flC)]) should be clarified. By defini- 
tion, 7 _1 (for the logarithmic landscape) and 2 - " 7 (for 
the linear landscape) are the transition probabilities over 
the barrier of the hierarchy level 7 separating two basins. 
To get the transition probabilities between the specific 
points x (located in one basin) and y (located in the 
second basin) — just these probabilities enter in the ki- 
netic equation on ultramctric trees — one should divide 
the transition probability over the barrier by the number 
of states in the basin, which is in our case 2 7 . 

The value of a can now be found straightforwardly 
Rewrite the transition rate g 7 as 

g 7 = (I.6I8)- 7 = 2" Q7 (11) 

From (fTTj) one gets 

a = « 0.6942... (12) 

In 2 
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One sees from (|T0|) that in the large-time limit (t — > oo) 
the decay of the survival probability is exponentially fast 
on logarithmic landscapes. Since the survival probabil- 
ity consists of additive contributions from all possible di- 
rected paths on the 3-adic tree, only the linear land- 
scapes, i.e. the "principal" Fibonacci sequences, give 
the major contribution to the survival probability in the 
large-time limit (t — > oo), leading to the following alge- 
braic decay: 

P(t) ~ tr v = £- ln2 / lns s ~ £- In 2/ In 1.618 _ £-1-44 Qg\ 

where s g = 1 + r g = (V5+ l)/2 ~ 1.618. Thus, we have 
v w 1.44. 



V. CONCLUSION 

Despite our value v w 1.44 is much closer to numeri- 
cally obtained critical exponent v ~ 1.4 -r- 1.5 than many 
other values of v found in former analytic approaches (ex- 
cept the one found in [l(| and [lH fl2| ) , we are far from a 
naive thought that our simple consideration resolves the 
problem of analytic computation of the survival proba- 
bility in the nonlinear dynamical equation (JTJ in some 
region near the chaos border. We only have reformu- 
lated some particular problems of chaotic Hamiltonian 
dynamics in terms of Markov dynamics on boundaries of 
ultramctric trees with transition probabilities prescribed 
by internal dynamics of the system. Recall the two main 
ingredients of our consideration borrowed from Q: i) 
for K /• Kg the average local exit time, r„, from the 
metastable island of the hierarchy scale n is proportional 



to the number of periodic orbits, q n ; ii) for n — > oo the lo- 
cal exit time, t„, grows exponentially with the hierarchy 
scale n. 

Note that our description does not contradict with the 
value v = 3 found in scaling analysis [H just at the 
boundary K = K g . Namely our consideration "smears" 
the scaling r„ ~ s™ to some region below K g where the 
" secondary" Fibonacci sequences starting from the points 
O2, O3, O4 in the Fig |3] come into the play. 

Let us end up this letter by saying that the conjec- 
tured approach offers a possibility to rise some interesting 
questions concerning the internal structure of the stan- 
dard mapping ([1]). For example, it would be desirable to 
check numerically the existence of the "secondary" Fi- 
bonacci sequences starting from the points O2, O3, O4 in 
the FigJ3] In case of their presence, one would be intrigu- 
ing to think about the "phyllotaxis" [26|, [27j in chaotic 
dynamical systems. By this conjecture we would like to 
attract the attention of researchers working in nonlin- 
ear dynamical and chaotic systems to the language [22[ 
developed for the description of stochastic processes in 
ultrametric spaces. 
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